Energy deposition in hard dihadron triggered events in heavy-ion colhsions 
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The experimental observation of hadrons correlated back-to-back with a (senii-)hard trigger in 
heavy ion collisions has revealed a splitting of the away side correlation structure in a low to in- 
termediate transverse momentum {Pt) regime. This is consistent with the assumption that energy 
deposited by the away side parton into the bulk medium produced in the collision excites a sonic 
Shockwave (a Mach cone) which leads to away side correlation strength at large angles. A prediction 
following from assuming such a hydrodynamical origin of the correlation structure is that there is 
a sizeable elongation of the Shockwave in rapidity due to the longitudinal expansion of the bulk 
medium. Using a single hadron trigger, this cannot be observed due to the unconstrained rapidity 
of the away side parton. Using a dihadron trigger, the rapidity of the away side parton can be sub- 
stantially constrained and the longitudinal structure of the away side correlation becomes accessible. 
However, in such events several effects occur which change the correlation structure substantially: 
There is not only a sizeable contribution due to the fragmentation of the emerging away side par- 
ton, but also a systematic bias towards small energy deposition into the medium and hence a weak 
Shockwave. In this paper, both effects are addressed. 

PACS numbers: 25.75.-q,25.75.Gz 
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I. INTRODUCTION 

The experimental observation of hadrons correlated 
back-to-back with a hard or semi-hard trigger hadron in 
Au-Au collisions at 200 AGeV has revealed a splitting of 
the away side correlation peak in a semi-hard momentum 
regime between 1 and 2.5 GeV [ll, 0, H, |j| wich is absent 
in p-p collisions where two back-to-back peaks appear. 
This means that the main strength of the away side cor- 
relation in Au-Au collisions in this momentum region is 
not found in the direction of the away side parton but at 
a large angle with respect to it. This angle is found to 
remain constant if the trigger momentum is changed and 
also for a variety of associate hadron momenta in the 
semi-hard regime. This observation can be contrasted 
with back-to-back correlations at hard trigger and asso- 
ciate hadron momenta well above 4 GeV 5] which show a 
reappearance of back-to-back correlations as seen in p-p 
collisions, albeit suppressed. 

This pattern has given rise to the idea that while energy 
loss of a back-to-back parton pair is responsible for the 
suppression observed at high Pt, the measurements at 
intermediate associate hadron Pt show how this energy 
is redistributed into the medium and may in fact show 
the recoil of the medium in the form of a hydrodynami- 
cal Shockwave ^] . Phcnomenological comparisons of this 
scenario with the data using the same Monte-Carlo (MC) 
simulation for energy loss and energy redistribution in 
Shockwaves found agreement with both the high Pt cor- 
relation pattern [3, Q and the low Pt peak splitting Q . 
A comparison with the measured 3-particle correlations 
[lOJ has also been made in the same framework [ll| , how- 
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ever remains somewhat inconclusive as to prove or dis- 
prove the existence of Shockwaves as the chief mechanism 
for energy redistribution. However, in [12l | an important 
difference between sonic Shockwaves and other conical 
emission mechanisms has been pointed out, i.e. the lon- 
gitudinal elongation of the shock cone due to longitudinal 
flow which should result in a large extension of the corre- 
lation signal in rapidity for a hydrodynamical excitation 
of the medium. 

This elongation is obscured in single hadron triggered 
correlation measurements due to the fact that the ra- 
pidity of the away side parton is not determined by the 
rapidity of the trigger hadron and all possible rapidities 
of the away side parton have to be averaged. However, 
if the trigger is a sufficiently hard back-to-back hadron 
pair, then the rapidity position of the away side parton 
is very constrained and the elongation should be observ- 
able. Unfortunately, requiring a hard trigger hadron on 
the away side introduces a bias towards small energy de- 
position into the medium. In addition, an away side par- 
ton emerging from the medium does not only produce 
the leading away side hadron (which is part of the trig- 
ger) but also subleading hadrons building up correlation 
strength along the jet axis also at intermediate Pt, thus 
obscuring any large-angle signal of a Shockwave by filling 
in the dip between the Shockwave wings with a back-to- 
back peak. In this publication, we aim at a discussion of 
these effects. 



II. THE MODEL 

We simulate hard back-to-back hadron production in a 
Monte Carlo (MC) model. There are three important 
building blocks to this computation: 1) the primary hard 
parton production, 2) the propagation of the partons 



through the medium and 3) the hadronization of the par- 
tons. Only step 2) probes properties of the medium, and 
hence it is here that we must speciiy details of the evo- 
lution of the medium and of the parton-medium interac- 
tion. The model is described in great detail in [^, [l3| . 
here we will just provide a short overview. 



A. Primary parton production 

The production of two hard partons fc, Hn leading order 
(LO) perturbative Quantum Choromdynamics (pQCD) 
is described by 



^ AB^kl+X 

dp^dyidy2 



ij 



Xlft/Aixi, Q^)x2.fj/Bix2, Q^) 



da'^- 



dt 



(1) 

where A and B stand for the colliding objects (protons 
or nuclei) and 2/1(2) is the rapidity of parton k{l). The 
distribution function of a parton type i in A at a mo- 
mentum fraction xi and a factorization scale Q ~ px is 
fi/Aixi, Q^)- The distribution functions are different for 
the free protons 14, 15] and nucleons in nuclei jig. 1 17| . 
The fractional momenta of the colliding partons i, j are 



given by xi^2 = -^ (exp[±yi] + exp[±2/2]). 
Expressions for the pQCD subprocesses ^^—^ — (g^ t, u) as 
a function of the parton Mandelstam variables s^t and 
ii can be found e.g. in [l^|. By selecting pairs of k,l 
while summing over all allowed combinations of i, j, i.e. 
gg,gq,gq,qq,qq,7]q where q stands for any of the quark 
flavours u, d, s we find the relative strength of different 
combinations of outgoing partons as a function of pt- 
For the present investigation, we consider a dihadron trig- 
ger at midrapidity 2/1=2/2 = 0. By MC sampHng Eq. ([1]) 
we generate a back-to-back parton pair with given par- 
ton types and flavours at transverse momentum px- To 
account for various effects, including higher order pQCD 
radiation, transverse motion of partons in the nucleon 
(nuclear) wave function and effectively also the fact that 
hadronization is not a coUinear process, we fold into the 
distribution an intrinsic transverse momentum kx with 
a Gaussian distribution, thus creating a momentum im- 
balance between the two partons as pxi + Pt2 = ^t- 



B. Parton propagation through the medium 

The probability density P(a;o,2/o) for finding a hard ver- 
tex at the transverse position rg = (2^0,2/0) ^-^d impact 
parameter b is given by the product of the nuclear profile 
functions as 



P{xo,yo) 



TAiro + h/2)TAiro-h/2) 

TAA{h) 



(2) 



where the thickness function is given in terms of 
Woods-Saxon the nuclear density pA{r,z) as 7a (r) = 



/ dzpA (r, z) . Rotating the coordinate system such that 
the near side parton propagates into the (— x) direction, 
the path of a given parton through the medium ^(r) is 
determined by its primary vertex Tq and we can compute 
the energy loss probability P{AE)path for this path. We 
do this in a radiative energy loss picture [1^, [2^ by eval- 
uating the line integrals 

/•OO /"OO 

c^,(ro,(/.)= / d^iO and (gi) (rg ,</.)= / d^m 
Jo Jo 

(3) 
along the path where we assume the relation 

q{^)^ K-2- e^/" (0 (cosh p - sinh p cos a) (4) 

between the local transport coefficient q{£^) (specifying 
the quenching power of the medium), the energy den- 
sity e and the local flow rapidity p with angle a between 
flow and parton trajectory [21 Le and p are taken from 
medium evolution models [23, UM ^ discussed in [g] . 
LUc is the characteristic gluon frequency, setting the scale 
of the energy loss probability distribution, and (qL) is a 
measure of the path-length weighted by the local quench- 
ing power. We view the parameter _ftr as a tool to account 
for the uncertainty in the selection of as and possible 
non-perturbative effects increasing the quenching power 
of the medium (see discussion in [3]) and adjust it such 
that pionic Raa for central Au-Au collisions is described. 
Using the numerical results of [22l |. we obtain 
P{AE; ojc, R)path for ujc and R = 2ujl/{qL) for given jet 
production vertex and angle 4>. In the MC simulation, 
we first sample Eq. ^ to determine the vertex of origin. 
For a given choice of cf), we then propagate both partons 
through the medium evaluating Eqs. ^ and use the out- 
put to determine P{AE;uJc,R)path which we sample to 
determine the actual energy loss of both partons in the 
event. 



C. Hadronization 

Finally, we convert the simulated partons into hadrons, 
provided that a back-to-back pair emerges from the 
medium after energy loss. More precisely, in order to de- 
termine if there is a trigger hadron above a given thresh- 
old, given a parton k with momentum px, we need to 
sample Ai^'^{zi,pt), i.e. the probability distribution to 
find a hadron h from the parton k where h is the most en- 
ergetic hadron of the shower and carries the momentum 
Pt = Zi -PT- 

In previous works [7[,[8| we have approximated this by the 
normalized fragmentation function Dk^h{z, Pt), sam- 
pled with a lower cutoff Zmin which is adjusted to the 
reference d-Au data. This procedure can be justified by 
noting that only one hadron with z > 0.5 can be pro- 
duced in a shower, thus above z = 0.5 the Dk^h{z,PT) 
and Ai^'^ {zi , Pt) are (up to the scale evolution) iden- 
tical, and only in the region of low z where the frag- 
mentation function describes the production of multiple 
hadrons do they differ significantly. 



We improve on these results by extracting Ai{zi,pt) 
from the shower evolution code HERWIG 25]. The pro- 
cedure is described in detail in [ij]. Samphng Ai{zi,pt) 
for any parton which emerged with sufficient energy from 
the medium provides the energy of the two most ener- 
getic hadrons on both sides of the event. The harder of 
these two defines the near side. The hadron opposite to 
it is then the leading away side hadron. For the present 
investigation, we require both to be in given momentum 
windows to count a dihadron triggered event. We average 
the energy loss on near and away side parton over many 
such events to determine the average energy deposition 
into the medium. 

In order to compute the correlation strength associated 
with subleading fragmentation of a parton emerging from 
the medium we evaluate A2(zi,Z2,Pt) (also extracted 
from HERWIG), the conditional probability to find the 
second most energetic hadron at momentum fraction Z2 
given that the most energetic hadron was found with frac- 
tion zi. This contribution to the strength of the away 
side correlation is competing with the Shockwave signal. 
Our way of modelling hadronization corresponds to an 
expansion of the shower development in terms of a tower 
of conditional probability denities Aj^{zi, . . . , z„, /i) with 
the probability to produce n hadrons with momentum 
fractions zi, . . . Zn from a parton with momentum px be- 
ing Il'2^^Ai{zi, . . . ZijPr)- Taking the first two terms of 
this expansion is justified as long as we are interested in 
sufficiently hard correlations. However, in the following 
we also consider situations in which the near side trig- 
ger momentum is rather hard 0(10) GeV, the away side 
trigger momentum is likewise hard 0(5) GeV, but with a 
substantial gap between near and away side to allow for 
energy deposition in the medium, but observe fragmenta- 
tion yield associated with this trigger in a regime where 
hydrodynamics is valid, i.e. 0{l) GeV. Since the di- 
hadron trigger forces the parton to high momenta, multi- 
hadron production at the low associate scale is likely. 
Consequently, we have to include the next terms in the 
expansion. A detailed numerical treatment is very com- 
plicated, however we estimate the next two terms as 



A3{zi,Z2,Z3,Pt) ~ A2{zi + Z2,Z3,Pt)9{z2 - Z3) (5) 

and 



Ai{zi,Z2,Z3,Z4,pT) «A2(zi + Z2 + Z3,Z4,Pt) 
X e{z2 - Z3)e{z3 - Zi). 



(6) 



This procedure explicitly guarantees energy-momentum 
conservation and preserves the correct ordering in hadron 
momenta inside the jet. For the results quoted in the 
following, we have verified that the results converge and 
that A4 is only a correction, and that hence the inclusion 
of further terms does not alter the result substantially. 



III. RESULTS 

In Fig. [1] we show the away side energy deposition into 
the medium created in central Au-Au collisions at 200 
AGeV as a function of the trigger momenta on near and 
away side for two different medium evolution models, a 
hydrdynamical code [23[ and a parametrized evolution 
model [2J| . This is the energy available to excite a shock- 
wave. Note that according to the phenomenological anal- 
ysis [3, |Tl|, ll3 a large fraction / = 0.75 (but not all) of 
the available energy actually excites a Shockwave. 
The energy deposition is always largest when the gap 
between near side and away side trigger momentum is 
maximal. There is some dependence on what model for 
the medium evolution is assumed to be valid, however 
some general trends remain robust: The energy deposi- 
tion is roughly a third of the highest (near side) trigger 
energy. The additional variation with the away side Pt 
is about 50%. 

On the other hand, if no away side trigger is required, 
typically all of the energy of the away side parton is lost 
to the medium [s], |9| . Since the parton energy is on av- 
erage roughly a factor two more than the energy of the 
leading hadron, requiring a dihadron trigger reduces the 
signal strength of the Shockwave by about a factor six as 
compared to a single hadron triggered event. 
Let us now compare the strength of the Shockwave cor- 
relation signal with next-to-leading and higer order frag- 
mentation of the away side parton. For this compari- 
son, we consider the associate momentum range of 1-2.5 
GeV where the PHENIX collaboration has first seen in- 
dications for a Shockwave [ij. As explained in detail in 
[3, [nl, [141 , we cannot reliably compute the precise mag- 
nitude of the Shockwave per-trigger yield in a given mo- 
mentum window, especially as long as the trigger is in a 
semi-hard regime below 6 GeV, as the yield is not only 
dependent on assumptions about flow in the medium, 
but also recombination/coalescence processes [2a, [23, [2^ 
need to be addressed below this scale. However, let us 
boldly assume that the per-trigger yield in single-hadron 
triggered Shockwave events scales with the average trig- 
ger momentum and based on this assumption extrapolate 
from the PHENIX data with a trigger of 2.5-4 GeV to 
the two fragmentation-dominated trigger ranges of 6-8 
GeV and 10-12 GeV considered in this publication (note 
that there is good evidence from STAR data [4] that the 
rise of the yield is in fact substantially slower with trig- 
ger Ft)- With this maximal assumption, the per trigger 
yield given the PHENIX acceptance in the 1-2.5 GeV as- 
sociate momentum window for a 6-8 GeV trigger would 
be 0(2.5) and for a 10-12 GeV trigger 0(4), and, again 
on the level of a rough approximation, reduced down to 
0(0.4) and 0(0.7) in dihadron triggered events due to 
the bias on energy loss. 

On the other hand, the per-trigger yield into the 1- 
2.5 GeV associate momentum window due to subleading 
fragmentation of the away side parton can be computed 
in our hadronization scheme using the approximations for 
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FIG. 1: Average energy deposition on the away side for different dihadron trigger momentum ranges. The left panel shows 
6-8 GeV momentum for the hardest hadron, the right panel 10-12 GeV. The i— axis shows different bins for the away side 
hadron momentum, the y— axis the corresponding energy deposition for two different models of the medium evolution (see [g] 
for details). 





1 1 


1 1 1 1 


- 


• neai- side P^ 6-8 GeV 






. near side P.J, 10-12 GeV 


- 




- 


i 


* 


- 




i 


i 




- 


1 


1 1 1 1 


- 




4-5 


5-6 6-7 


6. 



2 
1.8 

U 

DO 1.6 

&1.4 

Ij 

1.2 
1 



away side P bin [GeV] 

FIG. 2: Estimated per-trigger yield into the 1-2.5 GeV mo- 
mentum bin due subleading hadron production in fragmenta- 
tion of the away side parton as a function of near side and 
away side Pt range. 



^3 and A4 described above. Tlie results are summarized 
in Fig. H 

As seen from the figure, the yield is chiefly determined 
by the highest momentum scale (i.e. the near side trig- 
ger momentum) which is natural, given that this sets 
the overall energy available for hadron production in the 
jet. As the away side momentum scale is increased, the 
associated yield decreases. This is not unexpected, as re- 
quiring a larger fraction of the parton momentum to end 
up in the leading hadron, less momentum is available for 
subleading hadrons. 

However, the most striking result is that the expected 
per-trigger yields are of order 0(1), i.e. they are in fact 
by about a factor two larger than the upper limit for 
the per-trigger yields caused by the medium recoil due 
to the Shockwave. This means that if dihadron triggers 



are used to study Shockwave production, the dominant 
signal at midrapidity where the away side trigger hadron 
is observed is not the shock cone, but rather hadrons 
produced in NL fragmentation processes of the trigger 
parton. The Shockwave must then be observed as a cor- 
rection to this signal. Most importantly, a splitting of 
the peak with a dip at zero degrees and strength at large 
angles is not expected under these conditions. 



IV. SUMMARY 

We have discussed the expected changes in the correla- 
tion pattern seen in a hydrodynamical momentum regime 
when one goes from single hadron triggered events to di- 
hadron triggered events. The main advantage of a di- 
hadron trigger is that the rapidity of the away side par- 
ton is tightly constrained, thus a study of the medium 
recoil on the away side as a function of rapidity becomes 
meaningful. However, there are two effects which compli- 
cate the observation of the medium recoil substantially. 
First, by requiring a hard away side hadron, there is a 
significant bias towards events in which little or no energy 
was deposited into the medium. This reduces the energy 
available to excite a Shockwave, and hence the strength 
of the correlation by at least a factor six. 
Furthermore, once a hard away side hadron is detected, 
it is almost unavoidable that subleading, softer hadrons 
are created within the shower. This contribution is rather 
strong at low momenta and competes with the bulk recoil 
of the medium. We estimated here that it is at the posi- 
tion of the away side parton about a factor two stronger 
than the medium recoil. 

However, it is possible to eliminate the latter contribu- 
tion due to its different shape in rapidity: While any 
Shockwave signal is expected to be elongated in rapidity 



due to longitudinal flow, the jet cone due to fragmenta- 
tion in vacuum would not be elongated at all. Thus, by 
observing associate hadron production displaced in rap- 
didity from a hard dihadron trigger, a (weak) Shockwave 
signal should become visible without any contamination 
from soft hadron production in the jet. 
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